Large quasiparticle thermal Hall conductivity in the superconductor Bai_xKxFe2As2. 
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Both the thermal conductivity k^x and thermal Hall conductivity K^y in Bai_xKxFe2As2 display 
large peak anomalies in the superconducting state. The anomalies imply that a large hole-like 
quasiparticle (qp) population exists below the critical temperature Tc, consistent with either a 
highly anisotropic gap or a gap with nodes. The curves of n^x measured to 35 T fit well to a 
vortex-scattering expression. The qp mean-free-path inferred reproduces the observed anomaly in 
Kxy. Both channels provide a consistent estimate of a large qp population below Tc. 



The discovery 0, i, i, i, Q of superconductivity in the 
iron pnictides has galvanized intense interest in this new 
class of superconductors. The highest critical tempera- 
ture Tc achieved to date is ~55 K in the (1111) family. 
Superconductivity has also been found in the 122 iron 
pnictides based on FcAs and FeSe [(Sj. As in the cuprates, 
one of the key issues is the gap symmetry. Nuclear mag- 
netic resonance (NMR) relaxation experiments [7| are 
nominally consistent with nodes in the gap. However, 
angle-resolved photoemission spectroscopy (ARPES) ex- 
periments on Bai_xKxFe2As2 [1, Q favor isotropic gaps 
on the hole F Fermi Surface (FS), but the question of gap 
anisotropy is far from decided, especially on the sheets 
surrounding M . 

If the gap parameter A(k) is isotropic on each FS 
sheet, the population of Bogoliubov quasiparticles de- 
creases sharply below Tc (k is a wavevector on the FS). 
By contrast, if the nodes exist in A(k) (or if |A(k)| is 
strongly anisotropic), the quasiparticle (qp) population 
decreases quite gradually. In superconductors with un- 
conventional gap symmetry, the thermal Hall conductiv- 
ity Hxy bas proved to be a powerful probe for quasipar- 
ticles (qps). Unlike the diagonal thermal conductivity 
Kxx which is the sum of the electronic term Kc and the 
phonon term Kph, the off-diagonal term is purely 
electronic. Together, Kxx and Kxy have been used to 
probe extensively the qp densityand their lifetime in the 
cuprate YBaaCugOy (YBCO) 0, 
fermion superconductor CeCoIns 



121 and the heavy 
14| . By contrast, 
Kxx decreases roughly linearly with {Tc — T) below Tn in 
the clean s-wave superconductors Pb,Hg and Sn [l^, 
(data on Kxy are unavailable). 

We report detailed measurements of Kxx {T, H) and 
Kxy{T,H) on single crystals of Bai_xKxFe2As2 in the 
geometry with the field H||z||c and — VT||x. The lon- 
gitudinal and transverse temperature gradients STx and 
STy, respectively, were measured using chromel-alumel 
thermocouples. The 2 crystals studied have dimensions 
2 X 1 X 0.1 mm'^. At 10 K, the resolution achieved is 



STy - 10 mK. 



Measurements of Kxy were made to a max- 



imum field of 14 Tesla using thermocouples. To extend 
measurements of Kxx below 6 K where the thermocouple 
sensitivity falls steeply, we used matched RuOa; micro- 
sensors which are especially sensitive below 4 K (these 
measurements were performed to iJ = 35 T). 



In the Boltzmann theory approach, Kg from qp excita- 
tions is given by [H, [l^ 
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where E]^ — ^ A(k)2 + is the qp energy with ek the 
normal-state energy. Here, v(k) is the qp group velocity, 
/o the Fermi-Dirac distribution and T(k) the transport 
relaxation time. A recent treatment of Kxy applied to 
YBCO is given by Durst et al [13]. 
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FIG. 1: The thermal conductivity Kxx of single-crystal 
Bai_xKxFe2As2 at _ff = (circles) and at 35 T (triangles). 
The solid curve is the in-plane resistivity p in zero H (Tc = 
37 K). 

Figure [1] plots the T dependence of K^^iT, H) in Sam- 
ple 1 with iJ = (circles) and with 77 = 35 T (triangles) , 
together with the in-plane resistivity p (solid curve) . The 
curve for Kxx{T,0) is closely similar in Sample 2. As 
shown, Kxx{T,0) is nearly T independent between Tc {— 
37 K) and 100 K. Below Tc, it rises to a broad maximum 
that peaks near ^Tc. In a 35-Tesla field, the anomaly is 
almost completely suppressed. Both the anomaly profile 
and its field suppression are similar to features seen in 
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YBCO and CeCoIns. In these unconventional supercon- 
ductors, also rises to a large peak near iT^, reflecting 
a large qp population and a greatly enhanced (zero-field) 
qp mean- free-path ^o- 
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FIG. 2: Curves of the thermal Hall conductivity f^xy VS. H in 
Bai_xKxFe2As2 at T = 250^35 K (Panel a) and from 35^7 
K (b). Above Tc (37 K), K^y is H linear up to 14 T, but as T 
decreases below Tc, curvature becomes increasingly apparent 
(Panel b). At all T, the Hall signal is hole-like. 
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FIG. 3: The T dependence of the weak-field slope 
lim£f_o Hxy/H (solid circles with error bars). As T falls below 
Tc, f^'xyo displays a steep increase to reach a peak at 18 K, 
reminiscent of the behavior in YBa2Cu307. For comparison, 
the quantity AoKc (with A — 0.033, see below) is plotted as 
the dashed curve, and the electronic term Kc, inferred from 
Eq. [2j is displayed as open circles. 



The curves of K^y vs. H at fixed T are displayed in Fig. 
[21 Both above and below Tc, the sign of Kxy is positive 
(hole- like). Above Tc, Kxy is strictly linear in H (up to 
14 T). As T decreases from 100 K to Tc, the slope of Kxy 
vs. H increases gradually until Tc, where it undergoes a 
sharp increase. 

In the normal state, Kxy originates from the Lorentz 
force acting on the charge carriers (the Hall effect is also 
hole-like). Below Tc, the scattering of qps from pinned 
vortex lines possesses a right- left asymmetry which leads 
to a large Kxy [lOl . [T^ . [isj . The asymmetry originates 
from the circulation of the supercurrent around the vor- 
tex core and the Volovik effect. The qp Hall current 
initially scales linearly with the vortex line density ny — 
\B\/(f)Q, where B is the flux density and (pQ the supercon- 
ducting flux quantum. At large B, however, vortex scat- 
tering also severely reduces the mean-free-path of the qps 
(see below), resulting in a negative curvature in Kxy{H). 
In Fig. [lb, the steady increase in curvature is clearly 
seen as T decreases from 20 to 7 K. 

Focussing on the weak-if regime, we see that the 
sharp change at Tc in the qp Hall conductivity Kxy 
becomes apparent when we plot the weak-field slope 
K'^yQ = \m\H~,o Kxy/ H (solid circles in Fig. [3]). As T 
decreases from 200 K to Tc, K'^y increases slowly by a 
factor of ~3. At Tc, K'^yQ exhibits a sharp break in slope 
followed by a much steeper rise to a peak that is ~5 times 
larger than its value at Tc (we discuss the curves of AkcU 
and Kc below). 

We next turn to the diagonal term Kxx {T, H) , which 






(b) - 




22.5 K 




/ 17.5 . 




S;-^;/^ / 15 




12.5 


A 


^V-- 2-5 




1.2 



5 10 

HoH (T) 



10 20 30 

(T) 



FIG. 4: Curves of n^x vs. H at fixed T measured with ther- 
mocouples to 14 T (Panel a) and with RuOi, sensors to 35 T 
(b). Panel (a) shows in expanded scale the singular H depen- 
dence at T below Tc (37 K). Panel (b) displays measurements 
extended to = 35 T and to T = 1.2 K. In both panels, fits 
to Eq. [2] are shown as thin curves. 



provides quantitative estimates of the electronic term Kc 
and £o independent of Kxy (Fig. [4]). In the normal state, 
over the interval 250 to 40 K, the field dependence of 
undetectable with our sensitivity. Just below Tc, 
a weak H dependence becomes apparent (curve at 35 
K Panel (a)). As T decreases, this rapidly evolves to a 
singular \B\ dependence that characterizes the scattering 



3 



of qps from pinned vortices in type-II superconductors in 
the clean limit {Iq ^ where ^ is the coherence length). 
In Panel (a), the measurements of k^x were made with 
thermocouples to = 14 T. The measurements were 
extended to lower T (1.2 K) and higher H (35 T) by the 
use of matched RuOx micro-sensors (Fig. |4|d). 

The data in both panels of Fig 0] fit very well to the 
vortex-scattering expression (shown as thin curves) [lOl . 

MM 
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that, above Tc, accounts for less than ^^^xx^ a fraction 
similar to that in YBCO jj,]. 

The fits also yield estimates of a = which 
we display in Fig. O). As T decreases from 35 to 1.2 
K, a rises to 3x the value it had at Tc- This reflects 
primarily the increase in £o. According to Cleary p^ . 
the cross-section atr — With the estimate ^ ~ 20 
A (from Hc2 ~85 T), we find that a at 1.2 K corre- 
sponds to £o — 1,200 A. This supports our assumption 
that Bai_xKxFe2As2 is in the clean limit £q ^ ^. 

Since is proportional to a nominal picture of 
the T dependence of the qp population may be obtained 
from the ratio Kg/io- This quantity is plotted as open 
triangles in Fig. [5)3. Initially, Kg/lo is nominally T- 
independent, but decreases linearly with T below 13 K. 
The initial plateau is also seen in Ke/^o measured in 90-K 
YBCO 

The fits of Kxx vs. H to Eq. [5] has allowed us to 
determine Kg, Kph, and a independent of Kxy To show 
that these estimates are consistent with the thermal Hall 
conductivity, we note that k^^^q £q should share the 
same T dependence as an^ ~ £q, in the semi-classical 
approximation. In Fig. [31 we have plotted AksOl with 
the scale-factor A = 0.033 (dashed curve) to compare it 



with the measured values of K'xyo- Within the uncertainty 
inherent in K^j,g, the T dependences may be seen to track 
each other quite well, especially between 12 and 35 K. 
Hence, both Kxx and Kxy provide evidence for a large 
hole-like qp population that persists down to T ~ \Tc. 
Moreover, the values of a ~ inferred from Kxx give a 
consistent description of the T dependences of both Kg 
and k' 



a;yO ' 



where a{T) = £oatr/4'Q with atr the vortex cross-section ^ 
presented to an incident qp. Equation [5] expresses the 
additivity of the zero-i? scattering rate (~ £q^) and the § 
scattering rate introduced by vortices £~^ = atrny- Be- 
cause Hc2 >80 T, the condensate amplitude is nearly 
unaffected by the applied H. Hence it is safe to assume ^8 
that Ke{T) is nearly independent of H, and the dominant 
contribution to the observed H dependence arises from 
vortex scattering. We have assumed that the phonon 
term Kph has negligible H dependence, consistent with 
I, where q is the average phonon wavevector for 

r<T,. 

From the fits to the curves in Fig. 21 we have deter- 
mined Ke{T), Hph{T) and a(T). These values are plotted 
in Figs, andOD, respectively. As shown in Panel (a), 
the phonon term Kph (open triangles) decreases monoton- 
ically as T decreases below Tc- By contrast, the zeio-H 
electronic term Kg rises to a broad maximum at 15 K 
and then falls. As evident, the monotonic profile of Kph 
implies that the peak in Kxx (in zero H) is entirely as- 
sociated with Kg. In turn, the large peak in Kg demands 
a large qp population that survives to ^ ^Tr- We note 




FIG. 5: (Panel a) Comparison of n^x with He and inferred 
from fits of curves in Fig. [3] to Eq. [21 Above Tc, Ke is 



less than 



Below Tc, however, increases rapidly to 



account for the entire anomaly in the observed k^x- Panel (b) 
displays a = £o(Jtr/4'o (solid circles) and Ke/a (open triangles) 
obtained from the fits. With atr ~ 26 A, we estimate that 
£0 — 1,200 A at 2 K. The quantity Kc/a is a measure of the 
qp population below Tc. 

We next discuss how disorder and inhomogeneities may 
affect these measurements. As these samples are first- 
generation crystals, disorder effects are a matter of ex- 
perimental concern. Are the observed heat-transport 
anomalies an artifact of disorder? We argue that this 
possibility may be effectively excluded. As Kxy ^g, a 
decrease in Iq produced by additional scattering from im- 
purities and lattice disorder strongly decreases the 
superconducting state. The trend has been confirmed re- 
peatedly in earlier superconductors. In optimally-doped 
YBCO, the anomaly in Kxx increases 3-fold in size when 
we compare the earliest twinned crystals with untwinned, 
high-purity crystals grown in BZO crucibles at Univ. 
British Columbia |12l |. Values of Kxy show an even 
larger increase. In CeCoIns, the Kxx anomaly doubles 
in size between early crystals and later high-purity crys- 
tals T^. Instead of being the origin of these heat- 
transport anomalies, disorder actually suppresses them 
rapidly. We therefore conclude that they reflect the in- 
trinsic qp properties of the pure system. 

The value of ~1,200 A at 1.2 K - comparable with £0 
observed in the earliest generations of YBCO crystals - 
is surprisingly long. It provides a measure of the length- 
scale of the electronic homogeneity in Bai_xKxFe2As2. 
The long £0 weakens significantly arguments advanced 
to explain, for e.g., the power-law decay of the NMR 
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relaxation by invoking strong disorder. This system lies 
in the opposite (clean) limit. 

The abrupt increase in (T) at Tc is similar to that ob- 
served in YBCO fiol. [T^. In both systems, the increase 
seems too abrupt to be explained by the reduction in 
phase space for qp scattering caused by the opening of a 
strongly anisotropic gap. We raise the question whether 
this reflects the abrupt change of electron-electron inter- 
action caused by the onset of long-range phase coherence 
at Tc- This issue, never resolved in cuprates, has again 
become important in the pnictides. 

The symmetry of the gap in the iron pnictides has 
been the subject of considerable debate. The NMR re- 
laxation rate l/{TiT) in LaFeAs(Oi_a;Fa;) follows a 
power-law relaxation rate extending over 3 decades in T. 
Because an isotropic A would predict a Hebel-Slichter 
peak followed by exponential decay, the power-law relax- 
ation rate is evidence for nodes. In models in which A(k) 
is positive on the F-FS and negative on Af-FS (but oth- 
erwise isotropic), absence of the Hebel-Slichter peak may 
be explained, but the power law is problematical [20| . 

ARPES measurements on Bai_xKxFe2As2 reveal 2 
hole-like FS sheets around the point F and 2 elec- 
tron sheets around M, broadly consistent with LDA 
(local density approximation) band-structure calcula- 
tions [2l'|. Measurements by ARPES of the supercon- 
ducting A(k) d, Q yield a nominally isotropic gap on 
the 2 hole FS sheets at F. If the gaps on the F-FS sheets 
are indeed isotropic, we are left with the possibility that 
the qps inferred from heat transport reside on the elec- 
tron sheets (the ARPES resolution is insufhcient to rule 
out gap variation on the electron FS). The possibility 



of strong amplitude variation in A(k) when the FS ap- 
proaches the region in which A(k) changes sign has been 
proposed p^ . 

However, a new ARPES experiment has shown that 
the FS in Bai_xKxFe2As2 may be quite different from 
LDA predictions. Instead of the simple electron sheet at 
M, Zabolotnyy et al. [1^ recently observed a "propeller" 
shaped complex comprised of 5 FS sheets arrayed around 
the point X/Y. The 4 blades of the propeller correspond 
to hole ellipses while the cone corresponds to a circular 
electron FS. We reason that the existence of strong gap 
anisotropy on the hole sheets at X/Y will accomodate 
the large qp population, as well as the positive sign of 

As more experimental probes are brought to bear on 
the superconducting state in this interesting correlated 
multi-band system, we expect the picture to evolve. 
The properties of the quasiparticles extracted from heat 
transport will be valuable for understanding the pairing 
mechanism. 
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